Scaling properties of proton-nucleus total reaction cross sections 
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We study the scaling properties of proton-nucleus total reaction cross sections for stable nuclei 
and propose an approximate expression in proportion to Z' 2 ' 3 <r*° tal + 7V 2 ' 3 er*° tal . Based on this 
expression, we can derive a relation that enables us to predict a total reaction cross section for any 
stable nucleus within 10% uncertainty at most, using the empirical value of the total reaction cross 
section of a given nucleus. 
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Recently, we studied total reaction cross sections of 
carbon isotopes (N = 6 — 16) incident on a proton for 
a wide energy region [l|. In that work, using our nu- 
merical results, we found an empirical formula, Eq. (13) 
of Ref. [1] (see the Appendix). In this Brief Report, 
as an extension of the empirical formulae, we study the 
scaling properties of proton-nucleus total reaction cross 
sections for stable nuclei of the whole mass-number re- 
gion. Our new expression is found to be in proportion 
to Z 2 / 3 er*° tal + iV 2 / 3 cr£° tal . Although this is merely a pa- 
rameterization without any microscopic foundation, once 
one is able to accept this expression, we can derive a new 
relation (jJJ) that enables us to predict a total reaction 
cross section for any stable nucleus using only the empir- 
ical values of the total reaction cross section of a given 
nucleus at a given energy within 10% uncertainty. 

We assume here the following expression for the total 
reaction cross sections of a proton incident on a nucleus 
with the mass number, A(= N + Z): 
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where Z(N) is the number of protons (neutrons) 
in the target nucleus, Cp° tal (c*° tal ) is the proton- 
proton (proton- neutron) total cross section, and C(E)(= 
(C(E,A))a) is an energy-dependent coefficient to be de- 
duced from data. (• • ■) a implies the average value of the 
values for different mass numbers at each energy. Note 
that Eq. (JTJ is merely a parameterization without any mi- 
croscopic foundation and that the factorization of C(E) 
is crucial for later discussion. This expression can pre- 
dict the energy dependence of total cross sections of lead 
(Z = 82, N = 126) with acceptable uncertainties, but, 
unfortunately, the isospin dependence is not necessarily 
adequate for describing neutron-rich unstable nuclei. 

In Eq. ([1]), C(E,A) in the second line is replaced by 
C(E), because C(E,A) depends weakly on A. Figure Q] 
shows the values of ttC(E, A) as functions of the mass 
number for all available nuclei at six different energies [2j . 
The solid lines in each panel represent the values of 
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FIG. 1: Values of nC(E, A) as functions of target mass num- 
ber, A, for each energy, E, indicated by open squares. For 
convenience, we plot irC(E, A) instead of C(E,A). The solid 
lines indicate the average values. The dashed lines indicate 
standard deviations around each average value. The experi- 
mental data are taken from Ref. [21 . 



ttC(E) obtained by averaging ttC(E, A) over the whole 
mass numbers for each energy. The dashed lines denote 
the standard deviations of irC(E). These results validate 
the replacement of C(E, A) with C(E) in Eq. ((TJ). 

In Fig. [2j we summarize the values of ttC(E) as a func- 
tion of energy. Those values are also listed in Table Q] with 



LU 

o 



1.4 

1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 




200 



400 600 

E (MeV) 



800 



FIG. 2: Values of coefficient, ttC(E), as functions of proton 
incident energy, E. The solid curve shows the behavior of 
Eq. ©. 
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To facilitate the numerical calculations of ctr using 
Eq. (JTJ) , it would be useful to prepare the analytical form 
for C(E) as a function of proton incident energy, E. We 
adopt, for example, 



ttC(E) = oi -a 2 exp(-a 3 £; a4 ) 
x cos[a 5 £ ae ], 



(3) 



where _E is the projectile kinetic energy in units 
of MeV. The constants are given by ai=0. 909824, 
a 2 =2.329619, a 3 =0.4345765, a 4 =0.2655287, a 5 =0.075, 
and 06=0.6262386. The behavior of Eq. (j3|) is repre- 
sented by the solid curve in Fig. [2J 

Using the preceding setups, in Fig. [3] we show the re- 
sults of the formula with the data of <7r as a function 
of mass number, A, at two different energies. The solid 
curves show the results of Eq. ((T|) . The values of the coef- 
ficient, ttC(E), that provide the solid curves are listed in 



TABLE I: Values of coefficient, tyC(E), defined in Eq. ©, 
and their standard deviations. The curves in Fig. [3] are drawn 
using these values. 
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FIG. 3: Comparison of the numerical results obtained using 
Eq.fTJ) with the empirical data of total reaction cross sections 
for proton-nucleus reactions as a function of mass number, A, 
for E = 99 and 860 MeV. The solid curves are the numeri- 
cal results. The values of nC(E) are listed in Table [I] The 
experimental data are taken from Ref. Q|. 



Table |TJ The agreement of the numerical results with the 
empirical data is fairly good. Although we do not show 
the results for E = 40 MeV, one should be careful in the 
application of this formula at E — 40 MeV because of 
the large deviation from the data for A > 100. 

Also, we show the results of the formula with the em- 
pirical data of or as a function of proton incident energies 
in Fig. 01 The solid curves represent the numerical results 
of Eqs. (P) and ^ for the total reaction cross section of 
protons incident on two target nuclei, such as C and Pb. 

To show the precision of our formula, Eq. ([T]), we esti- 
mate the fluctuations of the data around the numerical 
results at each energy. We introduce the fluctuation de- 
fined by the standard deviations of C(E) divided by its 
mean values, that is, AC(E)/C(E), and plot them as a 
function of energy, E, in Fig.[SJ One can see that, except 
for the case of 40 MeV, the uncertainties in our numeri- 
cal results are less than 10%. That is, the feature of each 
nucleus affects the magnitude of the total reaction cross 
section by 10% at most. 

Let us derive the "parameter-free" expression of <7r us- 
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FIG. 4: Comparison of the numerical results obtained using 
Eqs.(frj) and Q with the empirical data of the total reaction 
cross sections for proton-nucleus reactions as a function of en- 
ergy, E. The solid curves are the numerical results of Eqs. |TJ) 
and ([3jl . The experimental data are taken from Ref. [2j (solid 
circles) and Ref. [3[ (open squares). 



ing Eq. (QJ, which is the major topic of this article. In 
Eq. ([1]), we introduced the coefficient, C(E) as being in- 
dependent of A within 10% uncertainty. If we accept this 
fact, we can estimate a total reaction cross section of a 
nucleus A\{= Z\ + N\) at a given energy using a known 
value of the total reaction cross section of another nu- 
cleus A2{= Z-i + AT 2 ) at the same energy. The expression 
becomes 
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where A implies the uncertainty coming from that of 
C{E). Here, A is defined by A= (S(E)) E , and S(E)= 
2AC(E)/C(E). (6(E)) E implies the average of 6{E) val- 
ues at various energies. We obtain A ~ 0.1. No fitting 
parameter appears here, but, at most, 10% uncertainty 
should be taken into account. 

In Fig. [6l we show the numerical results of Eq. (|4|) with 
the empirical data as a function of mass number, A, at 
99 and 860 MeV. For the reference nucleus, which cor- 
responds to A2 in Eq. (0|, we here choose a nucleus of 
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FIG. 5: Uncertainties in numerical results of or of Eq. ([T]) 
as functions of proton incident energy, E. The solid line is a 
guide for the eye. 



a medium mass number, such as 56 Fe, 64 Zn, and 64 Cu. 
For numerical calculations of Eq. ^ for 40, 99, 180, 
550 and 860 MeV, we choose a nucleus and the reac- 
tion cross section at each energy as ( 64 Zn, 1092±23 mb), 
( 64 Cu, 835±23 mb), ( 56 Fe, 662±19 mb), ( 64 Cu, 667±67 
mb), ( 64 Cu, 777±17 mb) and ( 64 Cu, 728±17 mb), re- 
spectively 0| . Here we neglect A. 

The expression (j4]) indicates that, if we have the empir- 
ical value of the reaction cross section of a given nucleus, 
we can predict the reaction cross section of other nuclei. 
The numerical results are consistent with the empirical 
data of stable nuclei. A possible way of estimating an 
unknown reaction cross section of a given nucleus at a 
given energy is to use Eqs. (p} and ([3]) or to use Eq. ((H). 
We do not need any parameter for the latter. The fluc- 
tuations in or coming from the feature of each nucleus 
affects the magnitude by 10% at most. 

In summary, we have introduced an approximate but 
new expression for proton-nucleus total reaction cross 
sections in terms of the number of protons and neutrons, 
and the total cross sections of proton-proton and proton- 
neutron reactions. This expression is applicable in the 
energy range from 40 MeV to 1000 MeV. Also, we have 
derived a simple relation for predicting the total reaction 
cross sections of a proton with any nucleus within 10% 
uncertainty, using the empirical values of the total reac- 
tion cross section of a given nucleus at a given energy. 
A better result is obtained if a neighboring nucleus is 
adopted. 
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where R(C) = 0.96±0.05 and N > 7. The above relation 
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FIG. 6: Comparison of the numerical results of Eq. ((4|) with 
the empirical data of the total reaction cross sections for 
proton-nucleus reactions as functions of mass number, A. The 
solid curves represent the numerical results obtained using 
Eq. (Q. The experimental data are taken from Ref. @]. 



Appendix A: Volume-type form 

We found that, in Eq. (13) of Ref. [l|, for all the carbon 
isotopes, the following relation is satisfied over the entire 
energy range: 
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suggests the following volume-type formula: 
a n (Z, N, E) = C{E) [Za^\E) + Na% ta \E)} . (A2) 

We find that this formula is able to estimate the 
proton-nucleus cross section in the mass number range 
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FIG. 7: Prediction using Eq. (|A2[) for the proton-nucleus total 
reaction cross sections at 99 MeV. The nuclei with 4 < A < 
30 fit the results with C(E) = 0.3523198 (solid curve) and 
those with 30 < A < 60 fit the results with C{E) = 0.27075 
(dashed curve). 



around b < A < fo+30, where b is any number. The appli- 
cability is limited to a certain restricted region in A. Fig- 
ure[7]shows the numerical results of Eq. (|A2[) . The nuclei 
with 4 < A < 30 fit the results with C{E) = 0.3523198 
(solid curve) and those with 30 < A < 60 fit the results 
with C(E) = 0.27075 (dashed curve). 
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